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We discuss the gravitationally interacting system of a thick domain wall and a black hole. We 
numerically solve the scalar field equation in the Schwarzschild spacetime and show that there exist 
scalar field configurations representing thick domain walls intersecting the black hole. 
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O ■ I. INTRODUCTION 

o ■ 



Topological defects arise during spontaneous symmetry breaking associated with phase transitions, and cosmological 
evolution of them is considered to have played an important role in cosmology (see e.g., Topological defects 

produced in the early universe might give us some information on high energy phenomena which cannot be reached 
by accelerator experiments, and those produced at the late time phase transitions also have attracted attention 
00 • as a potential source of the cosmic structures. To study the topological defects is therefore crucially important in 

| cosmology and elementary particle physics. 
t-H . In general relativity, topological defects have interesting features. Topological defects are extended and relativistic 
objects due to their large tension. In particular, the gravitational field produced by an infinitely thin domain wall 
shows repulsive nature HQ. On the other hand, although there are many studies about the properties of thick 
domain walls in the flat and de Sitter spacetime backgrounds [Q, little is known about the existence of thick wall 
ly-j , configurations on inhomogeneous, strongly curved background such as a black hole spacetime. Thus it is intriguing to 
study the gravitational interaction between two extended relativistic objects: a topological defect and a black hole. 

In most studies of defect-black hole system, topological defects have been approximately treated as infinitely thin 
and non-gravitating objects whose dynamics is governed by Nambu-Goto action. Within this context, the scattering 
""q I problem of a Nambu-Goto string by a background black hole has been studied in detail J^-^. In the domain wall 
Q-i case, recently M. Christensen, V. P. Frolov and A. L. Larscn considered Nambu-Goto walls embedded in the 

Schwarzschild black hole spacetime and found the static axisymmetric solutions. They showed that there exist a 
family of infinitely thin walls which intersect the black hole event horizon. However, little attention has been given 
to the system of a thick defect interacting with a black hole though a defect as a topologically stable configuration of 
. £h ! a scalar field has a finite thickness. 

Now we shall consider the validity of thin-wall approximation in the system of a topological defect and a black hole. 
In such a system, there are two characteristic scales: the thickness w of the defect and the black hole radius R g . In 
the case of the system of an astrophysical black hole with the mass ~ M® and a defect formed during a GUT phase 
transition, the thickness of the defect is much smaller than the black hole radius and therefore thin- wall approximation 
would be valid. However, it is not so hard to consider the situation that the thickness of defects cannot be negligible 
as compared with the size of a small black hole. Over the last few decades, many people have studied the formation 
of small black holes called primordial black holes (PBHs) and their cosmological implications. For example, studying 
the contribution of PBHs to cosmic rays enables one to place limits on the spectrum of density fluctuations in the 



early universe (see e.g. 12-14Q. On the other hand, it has been discussed the possibility of thick defects and their 
roles in cosmology, e.g. as a source of large-scale structure in the universe Q, or as a candidate for some kind of dark 
matter |HJ. For example, it is thought that the typical mass of PBHs which evaporate at the present epoch is 10 15 g, 
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so Rg ~ 10 _13 cm. When one considers the topological defects formed during a phase transition at < lOOMeV, such 
defects become thicker than the size of PBHs and thin-wall approximation is no longer valid. 

In this paper, we investigate the gravitational interaction between a domain wall and a black hole taking the 
thickness of the wall into account. We deal with scalar fields in the Schwarzschild black hole spacctime with 4> 4 and 
sine-Gordon potentials, which have a discrete set of degenerate minima. We explicitly show that static axi-symmetric 
thick domain walls intersecting the black hole do exist by numerical investigation. We consider a non-gravitating 
domain wall for simplicity. This test wall assumption might be valid when the symmetry braking scale of the scalar 
field is much lower than Planck scale as will be shown later by dimensional analysis. 

This paper is organized as follows. In section 2, we derive the basic equation and discuss the boundary conditions 
which represent the situation we want to study. In section 3, we show the numerical result. We summarize our work 
in section 4. We also discuss the validity of the assumption that the effects of gravity of domain wall can be ignored 
near the black hole horizon. Throughout this paper, we use units such that c = ti = G = 1 unless otherwise stated. 



II. THE BASIC EQUATION AND THE BOUNDARY CONDITIONS 

We consider a static thick domain wall in a black hole spacetime. The domain wall is constructed by a scalar field 
with self-interaction in a given curved spacetime. In what follows, we neglect the self gravity of the scalar field, as 
will be justified later. As a background spacetime, we consider the Schwarzschild black hole 

g=-(l-^jdt 2 +(l-^j cLR 2 + i? 2 (dtf 2 +sin 2 tfd^ 2 ). (1) 

For our purpose, we find that it is more convenient to work in the isotropic coordinates {t, r, 1?, ip\, where the new 
radial coordinate r is defined by 

M x 2 



R = r(l + -) . (2) 

We mainly consider the region outside the event horizon in this paper, which corresponds to r > M/2. In this 
coordinate system, the metric has a spatially conformally flat form 

r 2 + r 2 (dtf 2 +sm 2 i3d(p 2 )] . (3) 
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Let us consider a real scalar field <p with a potential V[4>], of which Lagrangian is given by 

C = -(-detg) 1 / 2 Qv0-V0 + y[0]) . (4) 

The equation of motion for cf> is 

V 2 <t> - dV/d<p = 0. (5) 

In this paper, we consider following two familiar types of potentials which have a discrete set of degenerate minima; 
the cj) 4 potential 

VM = - A i<f-ri 2 ?, (6) 

and the sine-Gordon potential 

VM = At, 4 [1 + . (7) 

Since the Schwarzschild spacetime is asymptotically flat, the asymptotic boundary condition for Eq. (|5|) would be 
given by the solution in the flat spacetime. The relevant solutions in the flat spacetime g = —dt 2 + dx 2 + dy 2 + dz 2 
are the static and plane-symmetric solutions 

0i (z) = r) tanh y/\/2r]z (8) 
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and 



62(2) = 77 [4 arctanexp(-\/A?72) — 7r], 



(9) 



for the potentials V\ and V2, respectively. These solutions represent domain walls in the flat spacetime characterized 
by the thickness of the wall 



l/VAry. 



(10) 



In the Schwarzschild background, the solution compatible with the above asymptotic boundary condition would have 
a static and axi-symmetric form <j) = <f>(r, •&). Then, the explicit form of the equation of motion (^|) becomes 



2r 



2r + M 



4 r 



dr 2 



8r 



d 



1 



(4r 2 - M 2 ) dr r 2 



d-d 2 
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d_ 
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This equation can be parameterized by a single dimensionless parameter 

e = M/2w, 

by introducing dimensionless variables 

p = 2r/M, $(p,0) =0(r,i?)/»7. 
In terms of these variables, Eq. (Ill]) becomes 



1 



<9/0 2 



2jo d 



p 2 \dd 2 
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(11) 



(12) 



(13) 



(14) 



where the dimensionless potential t/[$] = y[0]/A^ 4 is defined. U has minima at <f> = ±1 for 4 potential and at 
$ = (2n + l)7r (n = 0, ±1, ±2, • • •) for sine-Gordon potential. The parameter e is just a ratio of the horizon radius to 
the wall thickness measured in the asymptotic region, namely if e is smaller (larger) than unity, then the wall is said 
to be thick (thin) as compared to the size of the black hole. 

We shall confine ourselves to the case that the core of wall is located at the equatorial plane {i3 = it/2} of the black 
hole. The solutions without this assumption will be discussed in the separated paper [ fl6|| . Accordingly, we impose 
the Dirichlet boundary condition at the equatorial plane 



$1 



0=tt/2 



0. 



(15) 



Now it is sufficient to consider the north region, namely the solution in the south region can be obtained via $(p, i?) = 
— $(p, 7r — {ir/2 < 1) < 7r} of the spacetime. The regularity of the scalar field at the symmetric axis is give by the 
Neumann boundary condition 



9$, 
^9 



0=0 



0. 



(16) 



On the other hand, the boundary condition at the event horizon {p = 1} is given by the Neumann boundary condition 

(17) 

As is shown in Appendix, the condition (|l7| ) is the consequence of a natural requirement that the energy density 
observed by a freely falling observer remains finite at the event horizon. In practice, the region of numerical integration 
should be finite, so that we need an asymptotic boundary condition at p = /9 max for p max ^> 1. Taking into account 
the flat background solutions (||) and (||) , we impose the Dirichlet boundary condition 



$iIp=p„ 



and 



= tanh(2 1 ^ 2 ep max cos 1?) 



4arctanexp(ep max cos , i9) — 7r, 



(18) 



(19) 



for the <f> and the sine-Gordon potentials, respectively. 

In the next section, we numerically integrate the field equation (14|) using relaxation method under these boundary 
conditions at the equatorial plane (|l5|), at the symmetric axis (|l6), at the event horizon ( |l7| ) and in the asymptotic 
region (|q), (p~9[) for both the <j) 4 and the sine-Gordon potentials. 
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III. NUMERICAL INTEGRATION 



The scalar field configurations &(x, z) satisfying (|l4|) and the boundary conditions are shown in Figs. where x 
and z are the Cartesian coordinates x — psm-ff, z — pcosi?. Here we show the results in typical two cases; the e = 1 
case in which the thickness of the kinks ([l8]) , ( [l9| ) at /9 max is comparable to the Schwarzschild radius (Fig. [I] for <fi 4 
and Fig. || for sine-Gordon), and the e — 0.1 case in which the thickness of the kinks fll8|), ( |T^ ) at p max is one order 
of magnitude larger than the Schwarzschild radius (Fig. || for (f> 4 and Fig. |] for sine-Gordon). In both the cases, we 
obtain a domain wall solution as a kink of the scalar field at the equatorial plane z = 0. Particularly in the case 
e = 0.1, the black hole is enveloped in the core region of the wall. 

We also show the energy density E of the scalar field given by 



T t * 1 
~ Xr, 4 ~~ 2e 2 \p + 1 



a$\ 2 1 /a$ x 1 

dp J P 2 \ di9 



U[$] (20) 



in Figs. |^-|^, corresponding to Figs. |l]-^|, respectively. In all the cases, one can see that the configuration actually 
has a wall-like structure, namely the energy density is localized around the equatorial plane with a certain thickness 
corresponding to e. 

We shall comment on the computational domain and the grid spacing taken in our calculation. In order that the 
asymptotic boundary conditions ( |l8| ) and (jl^) make sense, p max must be large enough. In the above calculation, we 
take the computational domain which is fifty times as large as the horizon radius (i.e. p max = 51), and we carry out 
the integration on 500 x 90 grid (the grid spacing in p- and direction is 0.1 x horizon-radius and 1°, respectively). 
Then we clip the region {\x\, \z\ < 20}, where the above results are insensitive to the value of p max and the number of 
grid points. In fact, the results do not change when we extend the computational domain to p max = 101 with keeping 
the grid spacing, and the results differ from ones on the finer (1000 x 180) grid with keeping p max = 51 at most 1%. 
We also comment that the reliability of the numerical code is checked in the flat space case and it reprodeces the 
exact solutions (||), (^|) with accuracy of 10~ 2 . 

IV. SUMMARY AND DISCUSSIONS 

In order to answer the question whether or not scalar fields can actually form a topological defect in the vicinity 
of a black hole, we have numerically solved the equation of motion for real scalar fields with cj) 4 and sine-Gordon 
potentials, which have a discrete set of degenerate minima, in the Schwarzschild black hole background. In both 
4> and sine-Gordon potential cases, we showed that there exist the static axi-symmetric field configurations which 
represent thick domain walls intersecting the black hole. In particular, we studied the specific case that the wall's core 
is located at the equatorial plane of the Schwarzschild spacetime. We introduced the parameter e which characterizes 
the domain wall thickness compared to the black hole horizon radius; the smaller than unity e is, the larger the wall 
width is. We showed the domain wall solutions and their energy densities for e — 1 and e — 0.1 cases. In summary, 
we can say that a black hole is not an obstacle for scalar fields to form a domain wall configuration intersecting the 
black hole. 

One might wonder about our present results; one naively expects that the scalar field could have no static distribution 
around a black hole and inevitably fall into the horizon as usual objects do. However, a domain wall is a relativistic 
object with a large negative pressure (or large tension) whose magnitude is comparable to that of energy density. 
Furthermore in our study we examined the domain walls which are extended infinitely in space. Then we can 
understand that the domain wall is suspended from the asymptotic region and supported against falling into the 
black hole by its tension, so that the static configuration is realized. 

In our analysis, we assumed that the gravitational effect of the domain wall is negligible compared to that of the 
Schwarzschild black hole. We shall comment on the validity of this assumption. The energy density of a domain wall 
is given by 

GT t l ~ GX V 4 ~ -1 (JL ) , (21) 
ur \mpi/ 

where mp\ is the Planck mass and G is the Newtonian constant. On the other hand, the curvature strength, or 
gravitational tidal force, of Schwarzschild spacetime is estimated as 



GM 



i? 3 i? 3 



e- 



(22) 
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at areal radius R. Then a ratio of the gravity which will be produced by the domain wall to the gravity of the 
background black hole is given by 

ws ^?~K£) (^) ■ (23) 

When lj becomes much smaller than unity, the gravity of the domain wall is negligible compared to that of the 
black hole, and our test wall assumption becomes valid. From Eq. ( |2Sj ) we have u> ~ e 2 (rj/mpi) 2 near the horizon 
(i? ~ 2GM). Therefore, when we consider domain walls with the symmetry breaking scale being much lower than 
the Planck scale (i.e. 77 <C JWpi), we have w <C 1 and consequently our result gives a good description of shapes of 
gravitating thick domain walls near the black hole. We can also see from Eq. (f23|) that, in the thick wall (w > 2GM) 
case, if rj <C wipi, the test wall assumption is still valid even at R ~ w. 

In the asymptotic region (R 3> w), one may expect that the gravity of the domain wall is no longer negligible 
and changes the asymptotic geometry drastically. Bonjour, Charmousis and Gregory have recently investigated the 
spacetime of a thick gravitating domain wall with local planar symmetry and reflection symmetry around the wall's 
core fl7|| . They showed that the domain wall spacetime becomes spatially compact and has a cosmological horizon 
as de Sitter spacetime does. This suggests that, when a black hole exists and wall's gravity is taken into account 
in the region far from the black hole, the whole spacetime has a cosmological horizon and an axi-symmetric domain 
wall intersects both the black hole and the cosmological horizons as an equatorial plane in a Schwarzschild-de Sitter 
spacetime. This motivates us to study further the interaction between thick domain walls and black holes in, for 
example, Schwarzschild-de Sitter background. 

The domain wall solutions obtained here are thought to represent a possible final configuration of a gravitational 
capturing of a domain wall by a black hole. At present, it is far reaching for us to investigate a fully dynamical 
process such as the scattering and capture of thick domain walls by black holes. However, to get some insights into 
the problem, it is worth investigating the existence of the static axi-symmetric solutions which represent thick domain 
walls located away from a black hole. 

A cosmic string is also an extended relativistic object with large tension and thought to play more important role 
in cosmology than a domain wall does. Study of a gravitationally interacting system of thick cosmic strings and black 
holes is an interesting problem as a generalization of the present analysis. 
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APPENDIX A: THE BOUNDARY CONDITION AT THE EVENT HORIZON 



The tangent of a freely falling observer parameterized by its proper time t is = (di/dr, dr/dr, 0,0), and we have 



2r-AA 2 /d<\ 2 / AA 4 /dr x2 



1 -ft-v — ( JTi ld?J + l 1 + 5?J Is J ■ (A1) 



The quantity 



_ u {2r-M\ fdt\ ... 

is a constant of the motion, where £ M = (1,0,0,0) is the static Killing field. The energy density observed by this 
observer is 

T^vPu" = -(« ■ V^) 2 - (\v<t> • V0 + V[<l>]) . (A3) 
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Since we consider the st atic configuration and the spatial part of the metric is non-singular in the isotropic coordinates, 
the second term of Eq. (A3) is always finite. We have 



u ■ Vcf) 



M 
~2~r~ 



2r + M 
2r - M 



1/2 



dr 



(A4) 



for the static axi-symmetric configuration </>(r, Thus the requirement that Eq. (A3) is finite at the horizon is 
reduced to d(j)/dr = at r — M/2, or equivalently Eq. (|T^). 
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FIG. 4. Numerical solution of Eq. (^) with the sine-Gordon potential for e = 0.1. This solution expresses a domain wall 
configuration whose thickness is one order of magnitude larger than the Schwarzschild radius. The black hole is enveloped in 
the core region of the wall. 
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FIG. 8. The energy density E(x, z) of sine-Gordon scalar field for e 
shown in Fig. M 



0.1. This is calculated from the numerical solution 
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